An asymmetry model for the highly viscous flow 
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The shear flow is modeled in terms of local structural rearrangements. Most of these rearrange- 
ments are strongly asymmetric, because the embedding matrix tends to be elastically adapted to the 
initial state and to have a marked elastic misfit with regard to the final state. As one approaches the 
Maxwell time, the asymmetry becomes time-dependent, thus enabling the system to leave the initial 
state. The model explains the Kohlrausch behavior at the main peak in terms of the interaction 
between different local structural rearrangements. 

PACS numbers: 64.70.Pf, 77.22.Gm 
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I. INTRODUCTION 

The highly viscous liquid is a theoretical no-man's- 
land, with many theories and m0( Jelaii2iMi5iSiLSi2ii2JiiS 
competing for the explanation of the flow process. The 
elastic models^ claim a connection between the fragility 
(the abnormally strong temperature dependence of the 
main relaxation peak frequency) and the strong temper- 
ature dependence of the infinite frequency shear modulus 
G. But this claim has not yet solidified into a detailed 
microscopic theory. 

Experimentally, there is growing evidenc o 13 : 14 ' 15 : 16 for 
an intimate connection between the flow process itself 
(the primary or a-peak of the relaxation spectrum) and 
the faster relaxations three to six decades away. These 
faster relaxations appear sometimes as a secondary or 
Johari- Goldstein peak, sometimes just as an excess wing 
of the primary peak. The experiments support the idea^ 
that the a-peak is a result of the interaction between the 
"primitive" relaxations on its fast side. But again, this 
impression has not yet grown into a detailed microscopic 
theory. 

The present paper tries to do a first step towards a 
more detailed theoretical understanding. The asymme- 
try model proposed here is based on another recent ex- 
perimental finding! 7 - of a pronounced asymmetry of the 
secondary relaxation. These seem to be relaxation jumps 
from a lower energy minimum to a higher one, with an 
energy difference of about 4 fcsT g at the glass tempera- 
ture Tg. 

The asymmetry model ascribes the asymmetry to a 
structural rearrangement of an inner core, consisting of 
N atoms or molecules, with an N of the order of twenty 
to fifty. The structural rearrangement leads to a change 
of shape, which in turn leads to an elastic distortion, both 
of the core and of the surrounding matrix. The concept 
is a good starting point for a quantitative description of 
the flow process, not only because it provides a quantita- 
tive relation between the local flow event and its coupling 
constant to an external shear, but also because it allows 
for a simple treatment of the interaction between differ- 
ent relaxing entities. 

The following section II is divided in three parts. In the 



first, the asymmetry is explained in terms of the elastic 
distortion accompanying a structural rearrangement. Its 
average value of 4&bT is derived from the statistics of the 
structures. The second presents a proof of the 1/3-rule 
for interacting relaxation centers, adapted to the prob- 
lem. The third derives the shear response for a given 
energy landscape. Section III gives a comparison to ex- 
periment in two cases. Section IV contains the discussion 
and the conclusions. 



II. THE MODEL 

A. The basic picture 

The physical reason for a strong asymmetry of any 
structural rearrangement is rather obvious. Consider a 
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FIG. 1: (a) A structural rearrangement of three fiat molecules. 
Note the change of shape, (b) Potential energy, assuming an 
embedding elastic matrix adapted to the configuration of the 
left side. The dashed line shows the elastic contribution. 
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structural rearrangement of N atoms or molecules with 
an atomic or molecular volume v. Fig. 1 shows as an 
example the rearrangement of three flat molecules, in- 
volving a rotation of one of them by 90 degrees. The 
figure demonstrates one of the possible relaxation mech- 
anisms for rigid molecules, thus providing a physical ex- 
planation for the results of Johari and Goldstein^. Let 
the surrounding elastic matrix be adapted to the initial 
state, with no elastic strain of the volume. A transition of 
the volume Nv to a different stable structure changes the 
shape of the volume and leads to an elastic misfit, giving 
rise to an energy difference A between initial and final 
state. In the following, this mechanism is assumed to 
dominate the asymmetry of the double- well potentials to 
such an extent that one can neglect all other influences. 
The validity of this basic assumption will be discussed in 
Section IV. 

To put the concept on a quantitative basis, assume 
that the shape difference of the inner volume Nv is ei- 
ther a shear a (the shear angle between the two stable 
structures of the inner volume in radian) or a volume 
change Sv or a mixture of both. For a shear a and an 
infinitely hard surrounding medium, the asymmetry (the 
energy difference between final and initial state) A would 
be GNvef/2, where G is the (infinite frequency) shear 
modulus. But we consider a surrounding medium with 
the same elastic constants as the inner part. In this case, 
the medium distorts as well and one gets a lowering of 
the elastic energy by about a factor of two, slightly de- 
pending on the Poisson ratio of the substance^. For our 
purpose, the difference between the exact factor and two 
is negligible, so 



A = 



GNvei 



(1) 



and the coupling to an external shear e in the direction 
of e.j is given by 



dA 

de 



GNve, 



(2) 



One can bring the volume change NSv on the same 
footing as the shear by defining the scaled variable 



- 5 ,1 I 8B 
66 ~ VV3B + 4G' 



(3) 



where B is the infinite frequency bulk modulus and the 
factor results from the equilibrium between the com- 
pressed inner volume Nv, which cannot fully expand to 
its equilibrium value N(v + Sv), and the expanded sur- 
rounding clastic matrix. 

With this definition, the asymmetry A resulting from 
a structural transition of Nv with both shear and volume 
changes obeys again eq. {J) , with the square of given 
by a sum of the squares of its five shear components e\ 
to e§ and its volume change component e$ 



71=1 



(4) 



This means a given A defines the surface of a sphere in 
the six-dimensional distortion space. 

Let us next look at the statistics of the possible rear- 
rangements. Taking twenty to fifty atoms or molecules 
in the inner core, one can estimate the distortion needed 
for an asymmetry of 4kT. It is relatively small, less than 
a shear angle of ten degrees. So we are asking for the dis- 
tribution function of stable structures of the inner core at 
small distortions. It should be symmetric in the distor- 
tion, so it must start with a constant plus a second order 
term. This second order term will be negative for the 
volume change, because the average density is preferred. 
But there is no such preference for a shear. 

In the case of non-spherical molecules, each structural 
rearrangement will be accompanied by a change of molec- 
ular orientations, which in turn influence the change of 
shape. But the distribution function of accessible stable 
states must be still symmetric in the distortion, starting 
with a constant at small distortions. 

Therefore it is plausible to assume a constant proba- 
bility density of stable structures in the six coordinates, 
say p s jv in any of the five shear variables and Pvn i n 
the volume change cq. The number of states between 
ti and ei + dei is ir 3 pl N p V Ne^dei, because the surface 
of the six-dimensional sphere with radius a is 7r 3 ef. 
Since dA = GNveidei/2, one obtains a number of states 
p(A)dA between A and A + dA 



p(A) = 



n PsnPvn 
(GNv) 3 ' 



c N 



N 3 (Gv) 3 ' 



(5) 



where cn is a dimensionless and temperature- 
independent number, c/v will tend to increase with 
N, because a larger number of atoms or molecules offers 
more structural rearrangement possibilities 11 . N should 
be decidedly larger than six, because a recent boson peak 
investigation^ in terms of the Eshelby picture showed 
that the barrier in the case of six atoms or spherical 
molecules is still close to zero, the elastic restoring forces 
from the outside being of the same order as the ones 
from the inside. 

The partition function 



Z = 1 



cjv 



(GNv} 



r A 2 e fe s T dA 



(6) 



is easily evaluated. With the definition of the character- 
istic multi-minimum parameter /jy 



fN 



cn ( k B T 



N 3 V Gv 



the partition function reads 

Z = l + 2f N . 



(7) 



(8) 



The coupling 7 of a given structural rearrangement 
with asymmetry A to an external shear strain e 
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dA 
~de~ 



1 dA 

VGdei 



±y/AGNv/6, 



(9) 
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where the factor l/v6 stems from the average over the 
six possible directions of ej and the plus or minus sign 
depends on the orientation of e relative to e,, the plus 
sign occurring as often as the minus sign. Consequently, 
the term linear in e in the partition function vanishes. 
The second order term Z e results from the integral 



f N e 2 GNv 

where the function psq(A) is given by 

A 3 -a 



P«?(A) 



(k B T) 



-e k b j 



(10) 



(11) 



The A-dependence of the contributions to SG is pro- 
portional to A 3 exp(— A/k B T), a function with a maxi- 
mum at 3k B T and an average A of Ak B T. This is close 
to the asymmetry of 3.8 k B T found in the key experi- 
ment of Dyre and OlsenAI, a first important result of the 
model. 

To second order, the partition function with an applied 
strain e is 



Z = 1 + 21 



N 



f N e 2 GNv 
2k B T 



(12) 



With eq. (fl2j) . one can calculate the contribution 
5Gn of the volume Nv to the reduction of the shear 
modulus, using the free energy F — —k B ThiZ and 
SG N = (l/Nv)d 2 F/de 2 



SG 



N 



G 



1 + 2/ 



(13) 



From eq. (|13|) . it seems impossible to arrive at the 
breakdown of the shear modulus, 5G/G = — 1. However, 
as shown in the next subsection, the interaction between 
different relaxation centers increases their effectivity by 
a factor of three. Thus one needs only SG/G = —1/3, 
which is reached for /jv = 1. 

To understand the physical meaning of a characteris- 
tic multiminimum parameter 1, take a single asymmetric 
double well potential with A = Ak B T within the vol- 
ume Nv, with a coupling constant given by eq. f9|. In 
this situation f N « -5G/G = 1/6 cosh 2 2 w 0.012. This 
is the value where one crosses over from a single asym- 
metric double-well to a multiminimum situation. Thus 
the breakdown of the shear modulus requires hundreds of 
accessible elastically distorted structural minima around 
the undistorted ground state, a second important result. 



B. The 1/3-rule 

The 1/3-rule, first derived for asymmetric double-well 
potentials 2 ^, states that an ensemble of relaxation centers 
interacting via their coupling to the shear does only need 
an integrated strength of 1 /3 to bring the shear modulus 
down to zero. 



Here one needs a proof of the 1 / 3-rule for the multimin- 
imum case of n structural minima with energies U±, ..U n . 
Let us define 



Ui dui 



k B T de 



(14) 



where e is the external shear strain. The partition func- 
tion 



Z = y^exp(-Uj) 
l 

has the first derivative 

n 

Z' = -y^7»exp(-M t ) 
i 

and the second derivative 

n 

Z" = £ 7 2 cxp(- Uj ) 



(15) 



(16) 



(17) 



with respect to e. 

The second derivative of F — —k B T\nZ can be de- 
composed into pair contributions Fij with i =/= j 



„ k B T 
ij Z 2 



(7< - 7j) 2 exp(-it 4 - uj). (18) 



Next one considers the probability pt to find the system 
in the structural minimum i 



exp(— Ui) 

and its change under the applied shear strain e 

n 

. 7i - 7j 



z- 



■ exp(— Ui — Uj) 



(19) 



(20) 



On application of a small shear strain e at the time 
zero, the energy Uij transported into the heat bath by 
jumps between minima i and j in the equilibration of the 
n-minimum system consists of two terms, a term propor- 
tional to e and a second order term 



(2) 



k B T 
Z 2 



(a - l]) 2 exp(-u l - uj)e 2 



(21) 



The first order term has to disappear in the integral 
over the whole sample (otherwise it would not be sta- 
ble). The second order term is twice the free energy 
term F^e 2 /2. Since this holds for all pairs ij, it holds as 
well for the whole n-minimum system: The equilibration 
jumps transport twice the free energy change as heat into 
the heat bath. 

The proof of the 1 /3-rule proceeds by noting that the 
equilibrated n-minimum system is still a spanned en- 
tropic spring, able to deliver again the energy F"e 2 /2 
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to its surroundings. The equilibration of the whole liq- 
uid must convert also this stored energy into heat. Thus 
each relaxation center delivers three times this energy to 
the heat bath. 

Now consider an ensemble of volume elements Nv with 
Jn = lj which according to the preceding subsection 
brings the shear modulus down to 2/3 of its infinite fre- 
quency value, SG/G = —1/3. Apply again a small strain 
e at time zero. In the equilibration, the volumes trans- 
port two thirds of the total strain energy to the heat bath. 
The rest is lost in the full equilibration of the liquid. Since 
the energy transported into the heat bath cannot exceed 
the total strain energy, one concludes that the interaction 
between different relaxation centers increases their effect 
on the shear modulus by a factor of three. 

Note that the derivation of the 1 / 3- rule did not require 
the validity of the basic assumption of the preceding sub- 
section. 

The remaining task is to calculate the actual shear 
response for a given substance. This requires the quan- 
titative treatment of two nontrivial issues: (i) the dis- 
tribution of relaxation times, determined by the differ- 
ent barriers between structural minima and (ii) the time 
dependence of the asymmetry, reflecting the interaction 
between different relaxing entities. An approximation is 
derived in the next subsection. 



C. The shear response for a given energy landscape 

Let us consider the influence of the dynamics. In prac- 
tice, what matters is not the number of existing struc- 
tural alternatives, but rather the number of accessible 
structural alternatives. This depends on the relaxation 
times between the strain-free ground state and the sur- 
rounding strained excited states. If a given relaxation 
time, characterized 2 ^ by an energy barrier V, 



tv = t exp 



V 



(22) 



(with to = 1CP 13 s) allows one to see only one excited 
state, one is again in the situation characterized by a 
single double-well potential, even though there may be 
many other excited states with higher barriers around. If 
the temperature rises, or if one goes to longer relaxation 
times, they become visible. The condition fx = 1 is 
reached at a critical barrier V c . The corresponding decay 
time t c is the time at which the ground state transforms 
into an excited state of one of the surrounding states. 

To quantify the relaxation time influence, one defines a 
barrier density function f(V) in terms of the relaxational 
shear modulus decrease —5G/G = f(V)dV from all re- 
laxations with barriers between V and V + dV. f(V) is 
related to the conventiona l 2 ^ 23 relaxation spectrum H(t) 



The energy landscape of a given substance is charac- 
terized by the barrier density fa(V) of relaxation enti- 
ties, which reflects the barrier distribution between the 
ground state and the excited states. More precisely: If 
fNv(V)dV is the contribution to /jv from excited states 
between V and V + dV, then / (V) = /jvv(V)/3. For a 
given /o(V), V c can be determined from the 1/3-rule 

f (V)e X p(-Tv/T c )dV=^. (24) 

The final barrier density f(V) results from the interac- 
tion of the states fo(V), which leads to a time dependence 
of the asymmetry A. Let us consider the system in one 
of the excited states. As long as it does not jump back 
into the ground state, it is in a state of constant strain. If 
one holds a piece of matter under a constant shear strain, 
its elastic energy A decays as G(t)/G, where G(t) is the 
time-dependent shear modulus. 

The lowering of A increases the strength of a given 
relaxation, because the Boltzmann-factor exp(— A/fceT) 
increases. In order to calculate the enhancement for 
fo(V)dV, one has to do the integration of eq. (flU)). limit- 
ing oneself to barriers between V and V + dV. Therefore 
the enhancement factor is (A(0)/A(i)) 4 , as it would be 
for the full integral in eq. (jTDJ) 

h(V)G A 



f(V) 



(25) 



G{t v Y ' 

an equation which should at least be true at low V . One 
has 

I / f(V')dV. (26) 



G(t v ) 



G 







At low V, f(V) ~ fo(V); the enhancement factor in- 
creases with the integral of /(V) over V, which at the 
beginning is close to the integral of /o(V) over V. This 
suggests the following approximation 

fo(V)exp(-T V /T c ) 



f(V) 



where 



exp(— Ax — ax 2 ) 



f (V')ex P (-T v /r c )dV'. 



(27) 



(28) 



The term —Ax in the exponent of the enhancement factor 
of eq. (|27|) comes from the fourth power in eq. (j25|) , the 
coefficient a of the second order term is fixed by the 
normalization condition 



/ f(V)dV = l. 
Jo 



(29) 



H(r e 



V/k B T 



) = H(t v ) = Gk B Tf(V). 



(23) 



Trying different /o(V) to fit substances with and without 
secondary peak, one always finds a-values close to 6. 

The approximation provides an elegant way of de- 
scribing the crossover from the low-barrier region, where 
one sees the true energy landscape, to the interaction- 
dominated Kohlrausch behavior at the primary peak. 
The interaction is correctly taken into account to first 
order. 
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III. COMPARISON TO EXPERIMENT 



Eq. (|27|) allows to extract the energy landscape func- 
tion fo{V) from experimental data. This is easiest for 
substances with an excess win g 24 i 25 , which can be fitted 
with only three parameters: the infinite frequency shear 
modulus G, the decay rate r c and the slope 7 of the excess 
wing 



fo(V) = aexp 



fjV_ 
\k B T 



(30) 



where a is the normalization factor required by eq. (|24[) . 
Experimentally 2 ^, the exponent 7 was found to have val- 
ues around 0.2 in molecular glass formers. 

Having fo(V), one calculates f(V) from eq. (|2T)l . 
G'(uj) and G"{uj) are obtained from 2 ^ 



G 



G 



/WfV^Er (31) 



1 + w 2 r, 



< 32 > 



Fig. 2 (a) shows such a fit for triphenylethylene 2 ^ at 
260 K, 11 K above the glass temperature, with t c = 0.12 
s and 7 = 0.38. Fig. 2 (b) shows the best fit to glycerol 
shear data2£ about five degrees above the glass transition, 
with t c — 10.6 s and 7 = 0.14. 
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FIG. 2: Fit of experimental data (a) in triphenylethylene^ 
at 260 K with 7 = 0.38 (b) in glycerol at 192.5 K with 
7 = 0.14. In this second case, the fit also reproduces the 
excess wing of dielectric data— at 190 K. 



The fit in Fig. 2 (a) is close to the one of a Kohlrausch 
function with j3 = 0.57, but its wing has a slightly differ- 
ent slope, namely -0.38 instead of -0.57. This excess is 
even more pronounced in the glycerol case, where the 
a-peak itself is well fitted by a Kohlrausch law with 
(3 = 0.43. Instead, the fit with the asymmetry model pro- 
vides an excess wing with the much lower slope of 0.14, 
in good agreement with the one in the scaled dielectric 
data— at 190 K in Fig. 2 (b) (at high frequency, where 
one sees essentially fo(V), the dielectric e"(w) should 
show the same slope as G"{uS), as long as the coupling 
constants to dielectrics and shear maintain a constant 
ratio). 



IV. DISCUSSION AND CONCLUSIONS 

Let us first discuss the basic assumption of Section I. It 
assigns the energy difference between two stable inherent 
structures to the elastic distortion of the local structure. 
Obviously, this is not strictly true. One can imagine a 
structural rearrangement where the elastic distortion is 
zero. Then, the two energies might still differ, either 
because the energies of the two core structures differ or 
because their surface energy to the surrounding matrix 
differs^. 

However, neither the structural energy nor the surface 
energy supplies a quantitative explanation for the high 
average asymmetry found in experiment. Thus the basic 
assumption should be essentially correct, at least as an 
approximation; the other terms are not zero, but their 
influence on the asymmetry is probably small. 

Further support for the basic assumption comes from 
molecular dynamics results 29 which find a pronounced 
back-jump probability at shorter relaxation times, moti- 
vating a description in terms of metabasins, i.e. groups 
of inherent states which are only left for long relaxation 
times. Here, this is explained in terms of the elastic dis- 
tortion, which needs the Maxwell time to decay. 

A second question for discussion is the validity of the 
approximation for the shear response of a given energy 
landscape in the third part of Section II. It is a crude 
approximation, in particular for the region of the primary 
process. Nevertheless, it is a compromise which (i) works 
and (ii) illustrates an important fact: The shape of the cx- 
peak and in particular its slope on the high frequency side 
are only slightly influenced by the energy landscape and 
result mainly from the interaction between the relaxing 
units. One gets a Kohlrausch-/? of the order of 0.5 even 
if the underlying energy landscape has a much smaller 
(or even opposite) slope. Only at much higher frequency 
one starts again to see the true energy landscape. 

In this present version, one has to put in a secondary 
peak by hand, if there is one in experiment. This is not 
due to the crude way of treating the effects of the inter- 
action; whenever one has a secondary peak, it remains 
a peak in the glass phase, where the interaction effects 
disappear. 
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To conclude, the asymmetry model presented here de- 
scribes the flow process in terms of jumps from essen- 
tially undistorted ground states into elastically distorted 
states. The calculated asymmetry agrees with the exper- 
imental value. For high enough barriers, the return time 
becomes so long that the new state has time to relax into 
a new ground state. At this stage, the undistorted ground 



state can choose between hundreds of surrounding states. 
A recipe, based on a rather crude approximation, allows 
to calculate the shear response for a given energy land- 
scape. 

Valuable discussions with Roland Bohmer, Kia Ngai, 
Jeppe Dyre and Andreas Wischnewski are gratefully ac- 
knowledged. 
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